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Stability and surface diffusion at lithium-electrolyte interphases with 
connections to dendrite suppression 

Yalcin Ozhabes/'0Deniz Gunceler/'|^and T. A. Arias^ 
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This work presents an ab initio exploration of fundamental mechanisms with direct relevance to dendrite 
formation at lithium-electrolyte interfaces. Specifically, we explore surface diffusion barriers and solvated 
surface energies of typical solid-electrolyte interphase layers of lithium metal electrodes. Our results indicate 
that surface diffusion is an important mechanism for understanding the recently observed dendrite suppression 
from lithium-halide passivating layers, which were motivated by our previous work. Our results uncover 
possible mechanisms underlying a new pathway for mitigating dendridic electrodeposition of lithium on metal 
and thereby contribute to the ongoing efforts to develop stable lithium metal anodes for rechargeable battery 
systems. 
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I. INTRODUCTION 

Development of more efficient energy storage technolo¬ 
gies is needed to build a more sustainable future. Un¬ 
derstanding physical processes at the atomic scale on 
electrode-electrolyte interfaces is an important interme¬ 
diate step towards realizing this goal. Motivated by the 
desire to help enable many new applications, ranging 
from grid storage to long-ranged electric cars, researchers 
have also been using the tools of ab initio electronic 
structure to help develop better rechargeable lithium 

batteries 

The current state-of-the-art in rechargeable batteries is 
lithium-ion technology, where the presence of a graphitic 
anode host results in deadweight (carbon) to be car¬ 
ried along with the battery. Metallic anodes would be 
a better choice d ue to their increased energy density 
3860 mAhg~^)j ^^ l ^^ l but they suffer from localized 
nucleation while charging and form needle-like structures 
called dendritesW^^ Despite many years of concentrated 
effort, there are still many unanswered questions about 
the underlying physical mechanisms of lithium dendrite 
initiation and growth. This is partly due to the com¬ 
plex nature of the passivation layer, also called the solid- 
electrolyte interphase (SEI), that forms when the metal¬ 
lic electrode comes in contact with the elect rolyte.l^^^^^ 
Existing ideas and models on dendritic electrodeposition 
of lithium suggest that chemical inhomogeneities in the 
SEI layer result in spatially varying rates of deposition on 
the surface, which then lead to instabilities because any 
protrusion tends to concentrate electric field lines.l^^M^ 

Recent experiments have shown that the composition 
of the SEI layer has a dramatic effect on the performance 
of the battery cell. Based on our previous theoretical 
work,l^ Tu et. al. have managed to suppress dendrite 
growth in liquid and nanoporous electrolytes by passivat¬ 
ing the surface of the anode with lithium-halides.Like- 
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wise, researchers from Stanford have succeded in design¬ 
ing an interfacial layer from carbon nanospheres which 
improves cycling efhciency.l^ Understandably, these re¬ 
sults and many others stimulate a strong interest in 
studying the fundamental physical mechanisms within 
the SEI layer, and to date, many studies have investi¬ 
gated the bulk properties of these materials (e.g. bulk 
diffusion of Lithium) CnHUl 

The rich physics at the interface between anode and 
liquid electrolyte is much less understood, though there 
have b een very promising new developments in this field 
as well.Very recently, Jackie and GroB have pub¬ 
lished a comparative stuc^ of metallic Lithium, Sodium 
and Magnesium surfacesTheir DET calculations sug¬ 
gest that surface diffusion is significantly faster on mag¬ 
nesium metal than on lithium metal, which may be im¬ 
portant to understand why lithium forms dendrites while 
magnesium does not. While this work is very important, 
we believe (and the authors themselves also point out) 
that more investigation in this area is needed because 
that work does not address the presence and the effect of 
the electrolyte and, critically, the fact that metallic elec¬ 
trodes do not present pure surfaces to the electrolyte but 
rather complex non-metallic passivating layers known as 
SEE 

In this paper, we focus on the physical processes on 
the surfaces of various SEI materials for metallic lithium 
anodes. In particular, we provide an explanation of the 
physical mechanisms by which halogen additives (espe¬ 
cially E“) to the elect rolyte suppress dendrites and im¬ 
prove cycling efhcienc^ESHllESI utilize 

density functional theory to calculate surface cleavage 
energies and surface diffusion barriers for the most com¬ 
monly reported SEI materials in the liter at ure,SBHill ^^d 
then use these results to help understand the experimen¬ 
tally observed trends. Our hope is that such understand¬ 
ing will accelerate the development of novel anode ma¬ 
terials to improve batte ry perfo rmance. Our results sup¬ 
port the growing belieP^®^^!^ that anode materials with 
high surface energy and surface mobility are desirable. 
We also detail our previous claim,!^ recently supported by 
experiment,!^ that Lithium-halide SEI layers have these 
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FIG. 1. The percentage error in DFT lattice constants, plot¬ 
ted as a function of the experimental lattice constants. The 
solid black line represents exact agreement of theory and ex¬ 
periment. The dashed black lines represent ±2.5% deviations 
from experiment. 


desirable properties and that they are effective in sup¬ 
pressing dendrite growth on metallic lithium anodes. 


II. COMPUTATIONAL METHODS 

To perform first principles DFT calculations, we use 
the open source JDFTx softwar^^ which is based on the 
direct minimization of an analytically continued total en¬ 
ergy functional.!^ Ultra-soft pseudopotentialJ^ from the 
GBRV librar}!^ are generated using the Vanderbilt pseu¬ 
dopotential code.!^ To account for electronic exchange 
and correlation, we use the PBE flavour of generalized 
gradient approximation.!^ Throughout this methods sec¬ 
tion, we work in standard atomic units; i.e. bohr (B) 
for distances and hartree (H) for energies. All results in 
the sections below will be presented in more familiar SI 
units. 

For the Brillouin zone sampling of bulk units, we use 
a k-point grid of 4 x 4 x 4 which we determined by con¬ 
verging the total energy to a level of 0.1 mH. The energy 
cut-off for the plane wave basis was 20 Hartree which was 
also consistent with the same convergence threshold, a 
The nuclei were relaxed until the root mean square of the 
forces were below 0.1 mH/Bohr. 

To test our choice pseudopotentials as well as other 
calculation parameters, we calculated lattice constants 
of various Lithium SEI materials. The results, plotted 
in figure were satisfactory. (It is our belief that DET 
has the largest error in the lattice constant of LiOH be¬ 
cause LiOH has a layered structure where long-ranged 
dispersion interactions play an important role.) 

To calculate the surface formation energies, we create 
slabs varying in size from compound to compound. To 


determine the thickness of each slab, we take the number 
of layers that are necessary to converge the surface en¬ 
ergy to a level of 0.1 mH per unit cell. The vacuum layer 
(or solvent for fluid calculations) between the slabs is 20 
Bohr, allowing us to collapse the three dimensional k- 
point grid to a planar grid. The center layer is held fixed 
and the rest of the slab is relaxed with the same conver¬ 
gence criteria mentioned for the bulk calculations. Where 
needed, we use a truncated Coulomb kernel^ along the 
slab axis to prevent spurious electrostatic interactions be¬ 
tween slabs. Then we calculate surface formation ener¬ 
gies by taking the difference of the slab energy and the 
energy of equal number of formula units in the bulk and 
then dividing by the surface area. 

-^surf = ^^(-Uslab -^bulk) (1) 

To study surface diffusion, we consider the hopping 
process of an adatom that moves from one equilibrium 
position to the other and whose rate depends on the en¬ 
ergy barrier along that path. We determine the diffusion 
path by comparing the binding energies of the adatom at 
the high symmetry points along the surface, letting the 
adatom relax in the direction parallel to the slab nor¬ 
mal. Once we determine the endpoints of the diffusion 
pathway, we carry out a series of intermediate calcula¬ 
tions by putting the adatom in a series of sites equally 
spaced on the line connecting the two binding sites. Eor 
these calculations, the entire slab relaxes except the mid¬ 
dle layer, while the adatoms are restricted to stay on the 
plane perpendicular to the diffusion path. To minimize 
the interactions between periodic images of adatoms, we 
use 3x3 supercells. 

Standard plane-wave electronic structure methods 
have difficulty handling interfaces between battery elec¬ 
trodes and the electrolyte,!^ largely due to the need 
to thermodynamically sample the configuration space 
of the liquid electrolyte. As a result, there have been 
fewer ab inito investigations of lithium metal anode- 
electrolyte interfaces compared to, for example, the inves¬ 
tigations of t he b ulk properties of lithium intercalation 
compounds 

In principle, one can sample the configuration sp ace of 
the electrolyte using ab initio molecular dynamics 
and further accelerate the calculation using hybrid tech¬ 
niques such as QM/MM.!^ However, calculation of free 
energies with these methods are difficult and these ap¬ 
proaches often do not easily scale to the large number of 
materials we want to study. An alt ernativ e approach is 
that of continuum solvation models,!^^^^^^ where the in¬ 
dividual molecules in the liquid electrolyte are replaced 
with a continuum field, and thus free energies can be 
computed with a single density-functional calculation. 
Eortu nately, r ecent developments in continuum solvation 
model^^l^^H^ have made it feasible to efficiently study 
the solid-liquid interfaces for a large number of mate¬ 
rial/electrolyte combinations. 

In this work, we use a nonlinear polarizable continuum 
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FIG. 2. From left to right: rocksalt (Li-halides), anti-fluorite (Li20), zabuyelite (Li 2 C 03 ) and LiOH. Pictures were made 
using VESTAP 


model which models the electrolyte environment as a con¬ 
tinuous field of interacting dipolesThis approach can 
also capture dielectric saturation effects, which are im¬ 
portant near the surfaces of highly polar materials such 
as the Lithium SEI surfaces we consider here. In this 
nonlinear continuum model, the density profile of the 
electrolyte (s(f)) is computed self-consistently from the 
electron density of the surface slab (n(r)) as 


s{r) = erfc 


\n{n{f) / Tie) 
(tV2 


( 2 ) 


where a (= 0.6) determines the width of the transi¬ 
tion region that is set to be resolvable on the FFT grid 
and Tie is a (solvent-dependenfP^ critical electron density 
value that determines the location of the solute-solvent 
interface. For the non-electrostatic terms in the surface- 
electrolyte interaction, (such as cavitation entropy and 
long-ranged van der Waals) we make use of the effec¬ 
tive surface tension approximation,^^ which modifies the 
bulk (macroscopic) surface tension of the electrolyte to 
approximate these contributions. 


III. SURFACE ENERGIES AND DIFFUSION 

The formation of dendrites represents an increase in 
surface area. The thermodynamic perspective would thus 
indicate that SEI materials with greater surface energies 
should offer greater dendrite resistance. Enrthermore, 
dendrite nucleation may be driven by cracks in the SEI,^ 
so that a stable SEI with a high surface formation energy 
would offer resistance to this mechanism as well. 

However, the kinetics during electrodeposition might 
very well drive the system far from equilibrium. There¬ 
fore, one must also consider the mechanism by which 
surface energy would tend to suppress dendrites, namely 
surface diffusion. Eollowing the same train of thought, 
one expects that materials with fast surface diffusion, for 
example those with small diffusion barriers for adatoms, 
would be less likely to form dendrites. 

This section is broken into four parts: In the first part, 
we obtain order-of-magnitude estimates for those diffu¬ 
sion barrier heights which would make diffusion an im¬ 
portant process during electro-deposition. We find that. 


indeed, most SEI materials in the following sections do 
have diffusion barriers in the relevant range. In the sec¬ 
ond part, we investigate a binding-site switching mecha¬ 
nism that results in unusually low diffusion barriers for 
some lithium halides. In the third part, we summarize 
our results for the surface energies and diffusion barrier 
heights for many candidate lithium SEI materials, and 
discuss the trends we observe. We investigate multi¬ 
ple classes of materials: rocksalt halides (LiE, LiCl, Lil, 
LiBr), layered (LiOH) and multivalent (Li20, Li 2 C 03 ). 
In the final part of this section, we present an intriguing 
correlation between our diffusion barrier results and the 
experimentally observed time to short circuit. 


A. Comparing rates of deposition and diffusion 


cat: rates/main.tex: No such file or directory 
It is important to know whether, at experimental con¬ 
ditions, surface diffusion can be competitive with elec¬ 
trodeposition. We now probe this question with a simple 
order-of-magnitude analysis. We assume that diffusion is 
a random walk where the hop rate is given by the Arrhe¬ 
nius equation R = RQe~^, where kT is thermal energy, 
AE is the diffusion barrier and Rq is a base attempt rate 
related to the material’s phonon frequencies. To deter¬ 
mine the balance between surface diffusion and deposi¬ 
tion, let Qe be the charge of the arriving ions, J be the 
current density and a be the dimension of the surface unit 
cell. Under these conditions, we expect that dendrites 
form when the local current density becomes high enough 
to overcome surface diffusion, specifically when the mean 
distance traveled by a diffusing atom {d{t) = \fRta) in 
the (mean) time interval between two deposition events 
at the same site (t = q^j {Jo?)) is smaller than the grow¬ 
ing dendrite tip. In this simple scenario, the critical bar¬ 
rier at which the rate of electrodeposition equals the rate 


of diffusion, is given by AEc = 2/cTln 


J 


Experiments typically use current densities on the or¬ 
der of 1.0 mA/cm^, but owing to the highly non-uniform 
rate of deposition on the surface, especially near any 
protrusion, the relevant current densities can be orders 
of magnitude higher. Indeed, measurements of lithium 
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dendrite growtli^^^ have yielded the current densities 
at the dendrite tip in the 100 — 1000 mA/cm^ range. As 
for Rq, following previous studies,® we assume a standard 
phonon frequency of ^ 10 THz. For any microscopic pro¬ 
trusions in the order of 1-100 micrometers, we see that 
the critical diffusion barriers that result in competitive 
rates are of a few tenths of an electron volt. At 100 
mA/cm^ current density the critical barrier for a protru¬ 
sion of size 5 micrometers is approximately 0.1 eV, which 
would mean that for such a system diffusion on a surface 
with a barrier above this value would be too slow to miti¬ 
gate the growth of dendrites regardless of themodynamic 
unfavourability. In fact, we do find that surface diffusion 
barriers on typical SEI materials are in this range and 
so a detailed study of surface diffusion is necessary for 
understanding dendritic processes. Finally we note that, 
one must take this very simple analysis with care, as 
quantitative predictions other than order-of-magnitude 
estimates are likely to be beyond its capabilities. 


B. Change in the binding sites of Li-halides 

Below we find that of all the SEI materials we have 
studied, the lithium-halides are the most promising from 
the above point of view, with diffusion barriers ranging 
from 0.03 to 0.15 eV. The physical reason for these low 
barriers are as follows. Even though the bulk structure 
of all the lithium-halide materials we consider here is the 
same (rocksalt), the binding site for adatoms changes as 
the anion size increases. Eor halides with small anions (E 
and Cl) the binding site for the adatom is directly above 
the anion (“anion site”), and in the transition state for 
diffusion, the adatom sits between two anions and two 
cations (“in-between site”). See figure]^ (a) for an il¬ 
lustration. On the contrary, for large anions, the roles of 
these two sites are reversed and the binding site sits at the 
“in-between site” (slightly off-center) while the transition 
state has the adatom at the “anion site”. (See figure 
(b)) Now, the diffusion barrier is equal to the energy dif¬ 
ference between the binding site and the transition-site. 
Because of this, halide surfaces that are in the neighbour¬ 
hood of this change (LiCl and LiBr), where the sign of 
this difference changes, have very low diffusion barriers. 
On the other hand, LiE, which is far from this change, 
have a relatively larger barrier, but still small compared 
to non-halide SEI materials. 

The reason for the change in the binding site with an¬ 
ion size is steric interactions. Eor the two smaller halo¬ 
gens, the in-between site is too close to the two cations, 
which makes it energetically less favorable. However, as 
the anions get bigger, the distance from the in-between 
site to the cations increases and the adatom prefers to 
place itself in this halfway point where it can also inter¬ 
act with two anions simultaneously. Einally we note that 
even though the binding site changes as one goes from 
LiCl to LiBr, the diffusion path remains the same. 


C. Predictions for battery performance 


Eigure]^ summarizes all of our results, not only for the 
diffusion barriers for all of our surfaces (including halides, 
as well as the hydroxide, carbonate and oxide) but also 
for the corresponding surface energies. The data show 
that the presence of the electrolyte has significant impact, 
especially on the surface diffusion barriers. The surface 
energies of all ionic crystals go down (often by 5 — 15%), 
owing to the strong electrostatic interaction between the 
surface and the solvent. The surface diffusion barriers, 
however, change more dramatically, by up to as much as 
a factor of 2. The diffusion barrier for all materials but 
one (Li 2 C 03 ) decrease when the electrolyte is included 
in the calculation. These changes in diffusion barriers are 
significant because the rate depends exponentially on the 
barrier height and because these energy changes are on 
the order of several kT (which, at room temperature is 
approximately 0.025 eV). 

The data also suggest a strong positive correlation be¬ 
tween surface energies and surface diffusion barriers for 
most, but not all, SEI materials. The most severe trend 
breakers are Magnesium metal (black hexagon), lithium 
metal (black plus) and LiOH (blue square). Of these, 
the metals break the trend likely due to their very dif¬ 
ferent electronic structure, where Li is not in an oxidized 
(positively charged) state. LiOH likely breaks the trend 
because its layered structure and large intra-layer dis¬ 
tance cause it to have a very low surface energy along 
the z-axis. 

Li 2 C 03 , present in the SEI layer formed in the pres¬ 
ence of many commonly used electrolytes (propylene car¬ 
bonate, ethylene carbonate and others), is a less severe 
trend breaker. It has low surface energy and high dif¬ 
fusion barrier, both of which are undesirable quantities 
in a SEI material. Lithium halides, on the other hand, 
have lower surface diffusion barriers than Li 2 C 03 while 
also having either equal or higher surface energies. Our 
hypothesis, first put forth in an earlier work,® is that 
the low barriers may help explain the experimentally 
observed phenomenorP^ in which the formation of an 
lithium-halide SEI is effective in suppressing dendrites. 
We further hypothesize that these mechanisms may be 
relevant not only in experiments where the electrolyte 
has been seeded with a Li-halide crystal,!^ but also in ex¬ 
periments where other additives containing fluorine (e.g. 
hydrofloric acid or fluoroethylene carbonate]P^K^ have 
been used to improve stability and suppress dendridic 
growth. 

As for LiOH and Li20, which are also occasionally ob¬ 
served in experiments,!^ LiOH appears to be undesirable 
because of its low surface energy, whereas Li 20 appears 
to be undesirable because of its high diffusion barrier. 
However, the superiority of Li-halides over LiOH/Li20 
is not as conclusive because the halides are superior in 
only one of the two indicators. 

Einally, among halides, figure shows that the stabil- 
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FIG. 3. Surface binding energies versus its binding site. LiF is in the left, LiBr is on the right. The lowest energy binding sites 
are indicated with a green circle and the diffusion path is as shown with the green arrows. The black cross markers indicate 
our data points, the whole contour plot is generated using symmetries of the surface and cubic interpolation. 



FIG. 4. The surface diffusion barrier (x-axis) and the sur¬ 
face energy (y-axis) of various Lithium SEI materials. Hol¬ 
low data points represent surfaces in contact with vacuum 
whereas hlled data points represent the same surfaces in con¬ 
tact with the electrolyte (GH 3 GN, modeled with a polarizable 
continuum model). The black markers indicate the bare metal 
calculations by Jackie and GroB.^ 


ity (i.e. surface energy) decreases as one goes down the 
column of the peridic table, from F to Cl to Br to 1. This 
is likely due to some combination of the decrease in the 
electronegativity of the ions (which weakens the strength 
of the ionic bonds) and the steric interactions increasing 
the size of the lattice (which decreases the electrostatic 
stability of the lattice). Decreased stability is an unde¬ 
sirable property in battery materials as it tends to lower 


the voltage at which the surface breaks down. 

D. Comparisons with measurements of short-circuit time 


From a practical perspective, the ultimate figure of 
merit is the length of time before the battery short cir¬ 
cuits. Figure p lots this quantity, as measured by Lu, 
Tu, and Archer^ in a symmetric cell, as a function of our 
calculated diffusion barriers. Because diffusion is an acti¬ 
vated process, we choose a logarithmic scale for the verti¬ 
cal axis. Moreover, for Li 2 C 03 , Li20 and LiOH there are 
no separate data, instead, for this single measurement 
that we have, where there are no halide additives, the 
surface of the anode contains spatially-varying fractions 
of all three of these species. The observed breakdown is 
likely due to the most-dendrite prone among these three 
materials, and thus represents a lower bound for the life¬ 
time of each of the materials. Due to the notably lower 
barrier which we calculate for the hydroxide, we strongly 
suspect that the observed cycle lifetime is associated ei¬ 
ther with the carbonate or the oxide. Finally, because 
the barriers we calculate for these latter two materials 
are so close, we can not determine unambiguously which 
of these two materials fails first. Accordingly, the figure 
includes the experimental lifetime twice, once for each 
one of these two materials. 

The data in figurej^show a clear linear trend consistent 
with an Arrhenius-like behaviour. However, the slope 
does not appear to be exactly unity, demonstrating the 
importance of other factors beyond the diffusion barrier. 
Additionally, one point, lithium fluoride, appears to be 
an outlier, requiring additional explanation. LiF has the 
highest barrier among the halides but the second longest 
time to short circuit. We suspect that this deviation 
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Diffusion Barriers (in units of kT) 


FIG. 5. Battery short circuit times (from Lu, Tu, and 
Archer^^ plotted as a function of the surface diffusion bar¬ 
rier. The green line is the best fit line, excluding the one 
outlier. The blue line shows has the slope for an Arrhenius 
process. 


from the trend may be tied to it having the most stable 
surface among the halides (see figure |^, thus affording 
it additional stability. Lithium fluoride’s high surface 
energy derives not only from the high electronegativity 
of the fluoride atom, but also from its very small lattice 
constant and surface area per formula unit. (LiCl, the 
second smallest halide, has a surface area ^ 60% larger 
than LiF.) 

While it is not suprising that a single quantity (the 
diffusion barrier) is not enough to fully explain the ex¬ 
perimental lifetimes, we find it encouraging that there 
is such strong correlation between diffusion barriers and 
the experimental lifetimes. 


IV. CONCLUSION 

Recent experiments, prompted by earlier theoretical 
work,l^ confirm the success of Lithium- halide additives 
in suppressing dendrite growth,consistent with 
the findings of previous experiments with other fluoride- 
containing electrolyte additivesPrompted by this, 
we set out to explore more deeply the mechanisms of 
dendrite suppression at the atomic level. 

We performed density-functional calculations to deter¬ 
mine the surface energies and the surface diffusion barri¬ 
ers of lithium solid-electrolyte interphase materials. Our 
calculations show that a lithium-halide SEI layer results 
in increased stability of the surface (higher surface en¬ 
ergy) and faster diffusion along the surface (lower sur¬ 
face diffusion barrier for adatoms), both of which are 
likely important in explaining the above phenomenon. 
Furthermore, our results provide an explanation for the 
unusually low diffusion harries on halide surfaces, trac¬ 


ing this effect back to a change in the binding site. This 
change in binding site, in turn, is driven by the trends 
in the electronegativity and the sizes of the anions in 
the halide crystals. Finally, we observe a clear, approx¬ 
imately Arrhenius correlation between battery time to 
short circuit and surface diffusion barriers, that can be 
used as a guide to suggest new material systems to fur¬ 
ther enhance stability against dendrite formation. 

This work, which focused on solid-electrolyte inter¬ 
phase materials, leaves many directions yet to be ex¬ 
plored. In mitigating dendrite growth, multiple diffu¬ 
sion pathways are available, including surface diffusion, 
bulk diffusion, and diffusion at the SEI-metal interface. 
A study similar to this one, but for the SEI-metal inter¬ 
face, would be illuminating. Einally, with the potential 
importance of surface diffusion established, results from 
ab initio studies should be incorporated into macroscopic 
models of dendrite growth. 
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